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The vertical structure of T Tauri accretion discs 
IV. Self-irradiation of the disc in the FU Orionis outburst phase 
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Abstract. I investigate the self-irradiation of intensively accreting circumstellar discs (backwarmed discs). It is modelled using 
the two-layer disc approach by Lachaume et al. 1 2003 ) that includes heating by viscous dissipation and by an external source 
of radiation. The disc is made of a surface layer directly heated by the viscous luminosity of the central parts of the disc, and of 
an interior heated by viscosity as well as by reprocessed radiation from the surface. This model convincingly accounts for the 
infrared excess of some FU Orionis objects in the range 1-200 urn and supports the backwarmed disc hypothesis sometimes 
invoked to explain the mid- and far-infrared excesses whose origins are still under debate. Detailed simulation of the vertical 
radiative transfert in the presence of backwarming is still needed to corroborate these results and spectroscopically constrain 
the properties of intensively accreting discs. 
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1. Introduction 

It is now widely accepted that circumstellar discs accompany 
the process of star formation, all the more as such discs have al- 
ready been imaged in the millimetre (e x. DG Tau.lDutrev et alJ 
1 1994 and in the infrared (ex. HK Tau. Stand l eldt et al J 19981) . 
Matter accreting from the disc onto the star is supposed to build 
up solar-mass stars in the timescale of a Myr; yet recent stud- 
ies of star forming regions indicate way too low accretion rates 
among T Tau ri stars (TTS) —in the range 10" 10 -10" 7 M Q /yr, 
as shown by [G ullbring et al.l dl998l) in the Taurus region and 
iRobberto et alJ d2003l) in the Trapezium. Part of the matter is 
supposed to be accreted in earlier phases of high accretion but 
young stellar objects (YSOs) are then so embedded that their 
disc cannot be optically observed ("see lAndrlll994 for the 
stages of evolution). The other part might also be accreted dur- 
ing brief, periodic phases of intense accre tion of otherwise qui- 
escent TTS ( Hartmann & Kenvonlll996l) : FU Orionis objects 
(FUors), which feature accretion rates up to 10~ 4 M o /yr and 
which are seen to undergo an increase in luminosity o f mor e 
than 4 visual magnitudes in a few years faerbidl 196611 19771) . 
may represent such a phase. 

FUors are convincingly modelled by a self-heated vis- 
cous accretion disc (active disc) that overwhelms the stel- 
lar light and the properties of FUor outbursts are a lso well- 
studied (Hartmann & Kenvon 1985: iBell et alJll995l) . Yet ac- 



tive disc mod els fail to explain the mid- and far-infrared ex- 
cess 

(eg. IS 

imon &.Tovcel fl988) as well as weak silicate fea- 
tureat 10 [xm sometimes present in emission ( lHanner et alJ 
1998); these properties seem more typical of an irradiated disc. 
Altern ative models have tackled this problem: Lodat o~& Bertinl 
(12001 1) showed that the self-gravity of the disc can trigger an in- 
stability that produces additional warming and were able to re- 
produce the SED. However, this model cannot account for the 
silicate feature, which requires a temperature inversion at the 
surface of the disc. A more pleasing appr oach is the presenc e 
of a circumstellar envelope proposed by Ed ams et al. (IS 
on one (observational) hand, Kenvon & Hartmann (1991) suc- 
cessfully fit the spectra of two FUors with such an envelope. 
On the other (theoretical) hand, the envelope could serve as a 
reservoir of infalling material, replenishing the disc between 
outbursts. On the third hand, the envelope can account for a 
part of the extinction observed in these objects. 

iKenvon & Hartmannl Jl99ll) : EU {l999) also proposed 
that the inner hot parts of the disc are bright enough to heat 
up the outer ones (backwarming) and produce an irradiated 
disc-like SED. Their simulation uses a black-body disc model, 
ie. without vertical temperature profile, and convincingly ex- 
plains the order of the excess at 30-100 ^im, yet it fails to 
reproduce the SED at 1 0-30 um and does not predict a sili- 
cate feature in emission.lMalbet & Bertoutl(ll99ll) : ICalvet etal] 
dl99ll) : IChiang & Goldreichl i ll 9971) showed that irradiation can 
produce a hot disc surface radiating at shorter wavelengths 
and account for emission features ( CO bands, silicates). Using 
their results, Lacha ume et al.N2003l) developed a two-layer disc 
model in which irradiation by a central star and viscosity are 
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irradiation: 
X = 1 along bearru 




H(r) 




Fig. 1. Geometry of a backwarmed disc. The incidence angle 
<p of the incoming radiation onto the surface and the angle 6 at 
which the inner regions are seen from it depend on the profile 
of the surface H(r). 

taken into account. In this paper, I shall use this model with 
the disc itself as a source of radiation and investigate the back- 
warming in FUors. 

2. The model 

I use the tw o -layer approach described in detail by 
lLachaume et alJ J2003h : the surface of the disc is directly ir- 
radiated by the incoming radiation with an optical thickness 
of unity in the visible along the slanted path of the incident 
beam. Its vertical optical thickness in the infrared, to which 
radiation is mostly reprocessed, is then much smaller than one. 
Therefore, it is at a higher temperature than its effective temper- 
ature. The interior of the disc is both backwarmed by the sur- 
face and heated by vis cous dissipation, modelled with the a vis- 
cosity prescription by Shakura & Sunvaev ( 1973). A scheme of 
the model is presented in Fig.[T] 

2.1. Irradiation 

Since most of the flux originates from the first few solar radii 
of the disc, the irradiation source is seen as a central, infinitely 
flat ring from the the outer regions of the disc (1-100 AU). The 
radial flux distribut ion is, acc ording to the standard disc model 
bv lShakura & SunvaevMl973l) . 



nr 5 



(1) 



where M* is the mass of the star, ro < R m i n a constant, and /? m ; n 
the inner radius of the disc. The value of ro depends on the 
torque exerted on the inner rim of the disc. If the latter extends 
down to the star, the standard model assumes that the former 
is linked to the deceleration of the matter falling from the su- 
personic disc to the subsonic equator of the star, and ro « 7?*, 
the stellar radius. However, the torque exerted on a disc with an 
inner gap remains unknown; it should depend on the magnetic 
coupling between the star and the disc. Without clear theoret- 



ical or observational hints, I decided to assume ro 
the standard model. 

The total luminosity of the disc amounts to 



R-k, as in 



^disc — 



= (3 - 2m) u 



2nrF(r) dr 
2 GM*M 



(2) 



4R* 



where u = y/R*/R m i n < 1. At a given location on the disc's 
surface, the heating flux, ignoring the shadowing of some parts 
of the disc by the star, is: 



Fi n = (2/3 + cos 6) sin sin ( 



(3) 



r 2 +H 2> 

where 8 is the angle at which the heating ring is seen, <p the 
incidence angle onto the surface of the disc, and H the altitude 
of the surface above the mid-plane. In Eq.|3]the first term mod- 
els the "limb-darkening", while the second and third ones are 
projection effects. With an infinitely flat ring 

H 



= Arctan ■ 
r 

<b = Arctan - i 

\dr / 



(4) 
(5) 



The characteristic temperature of the irradiation is needed 
to determine the optical thickness of the surface layer — it de- 
pends on the penetration of the incoming radiation. We use 
the flux-weighted average of the viscous effective temperature 
7v(r): 



-—[' 

^disc JR., 



27T J F(r)7V(r) 4 dr 



45 - 80m + 36m 2 6 GM*M 



(6) 



3 - 2m 
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2.2. Location of the surface 

Chian g et al 1 (12001 1) provide a method to consistently find the 
location of the irradiated surface in an irradiated two-layer disc. 
However, the method is more difficult to implement in models 
that also include viscous heating, in which viscosity-dominated 
regions may present self -shadowing. Instead, we use d a simpler 
approach proposed by Chiang & Goldreich and state 

that irradiation-dominated regions meet 



Hcc h 



(7) 



where h is the vertical pressure scale-height of the disc. It ap- 
pears that such a prescription results in a surface profile H(r) 
close to that predicted by th e self-consistent m ethod in the 
irradiation-dominated regions ( lLachaume et all2003l) . Though 
one expects it to fail in viscosity-dominated regions, the impact 
on the disc structure is low, since these regions are by definition 
not affected by the irradiation. 

We now look for a self-similar solution H/r = r 7 . Within 
the approximation H(r)/r «c 1, Eq.[3]reads 

^/3(?)(<r-.>?)i 



oc h 2 r 4 . 



(8) 
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Table 1. Fundamental parameters of FUor models: accretion 
rate, inner radius of the disc, flaring H/r, type of dust opacity 
Kx, and extinction A v . 





M t M 

(M|/yr) 


^min 

(Re) 


H/r 
(at 10 AU) 


K X A v 
(ref.) 


FU Ori 


3.2 x 1(T 5 


4 


0.20 


i 1.3 


VI 057 Cyg 


2.0 x 10" 5 


2 


0.35 


i 4.5 


V1515 Cyg 


1.3 x 10 -5 


5 


0.35 


h 2.8 


ZCMaB 


1.4 x 10" 5 


2 


0.20* 


i 1.8 


i: inhomogeneous aggregates (Hennina & Stoanienko 1996) 
h: homogeneous aggregates (Hennina & Stoanienko 1996) 
* : not constrained 



The pressure scale height is obtained from the vertical gravity 
field of the star: 

^"a/SP- (9) 

where 

T c = (F irr /2) 1/4 (10) 

is the temperature of the interior of the disc in irradiation- 
dominated regions. The last two equations, together with H oc 
h, yield 

Hocr 3l2 F 1/& . (11) 

in v ' 

By subsituting Eq.[8]into Eq. I1 ll one finally derives 

H/rccr 113 . (12) 

This determination does not provide a estimate for the pro- 
portionality constant, which I leave as a free parameter in the 
model. 

3. Results & Discussion 

3.1. Strategy 

I apply this model to unembedded FUors, since more embed- 
ded objects cannot be properly described without considering 
reprocessing of the disc emission by the surrounded dusty ma- 
terial, which is beyond the scope of this paper. As a crite- 
rion for "unembeddedness" I choose a decreasing SED in the 
infrared, which leads me to consider three confirmed FUors, 
FU Ori, VI 057 Cyg, and V1515 Cyg, an d one FUor can- 
didate, Z CMa iKenvon & Hartmannlll99ll) . Their unembed- 
dedn ess is confirmed by their lim ited extinction Ay < 4 mag 
ICohen & Kuhilll979t bell et al.lll995l) Photometry data are 
taken from IKenvon & Hartmannl ill 99 lb for VI 057 Cyg and 
V1 515 Cyg, fromlThiebaut et alJ ll 995b for Z CMa, and from 
the lGezari et alJ dl999h catalogue for FU Ori. 

The models have been computed using wave- 
length-dependent grain opacities for dust aggregates by 
Henning & Stoanienko ( 1996), needed for the determination 
of the optical thickness of the surface layer and for the com- 
putation of the SEP. A t higher temperatures, I used grey gas 
opacities bv lBell & Lin ( 1994). Model parameters R min , M+M, 
and Ay can be determined by first fitting an active disc model 
to the SED in the range 0.3-3 um; they do not depend on the 
strength of the mid-IR excess to be modelled. The constant 



H/r (10 AU) is then adjusted to fit the backwarmed disc 
model to the mid- and far-IR excess. M+M and H/r (10 AU) 
are unambiguouisly constrained (by the 1-3 um flux and the 
10-50 um SED respectively), while R m i n and Ay are derived 
from the SED at A < 0.5 um. 

In all models, the viscosity parameter is a = 0.1 and the 
albedo of the surface J[ = 0.2. a has almost no influence on 
the SED at A < 100 urn because these discs are optically thick, 
and the presence of the albedo increases the visible and near- 
IR SED by a few percent while decreasing the mid- and far- 
IR excess by up to 20 %. Z CMa and FU Ori, which exhibit 
observational hints of not beeing seen pole-on, are assumed to 
be seen with i = 40°, other targets are modelled as pole-on 
(i = 0°). The inclination has little influence for angles i < 45° 
and mostly scales down the visible and near-IR SED by a factor 
cos f. 

3.2. Results 

The fits are displayed in Fig. |2] and their fundamental param- 
eters in Table For comparison's sake a standard viscous 
disc model and a black-bod y backwarmed disc, that is with 
no superheated surface, as in Kenvon & Hartmann 11991b . are 
shown along with the two-layer model fit. It appears that the 
black-body model gives a slightly higher far-infrared excess 
and a smaller mid-infrared excess than the present model and 
does not produce the 10 um silicate feature. 

3.2.1. FU Ori 

Up to 10 (im the SED of FU Ori can be modelled by a stan- 
dard viscous disc (Fig.|2] dotted line of the top-left panel). The 
accretion rate MM t = 3.2 x 10~ 5 Mg/yr corresponds to typi- 
cal values in the literature (« 3/c osiM^/yr where i is the in- 
clination of the dis c according to Hartm ann & Kenvonlll996l 
Mal bet et af]ll998l) . Nevertheless, the 5-12 um SED is about 
30% weaker as expected from the fit of 0.5-3 urn fluxes. 

A black-body backwarmed disc (see IKenvon & Hartmannl 
Il99ll) can explain the far-IR excess but it either overestimates 
the 30-100 um fluxes (same figure, dashed line) or underes- 
timates the 20-30 and 100-200 urn SED (not shown). If the 
backwarmed disc model comprises a hot surface resulting from 
direct irradiation, then the SED can be convingly fitted in the 
range 20-200 um (solid line) with H/r * 0.2 at 10 AU from 
the central star. This fit, however, overestimates the M and N 
fluxes by a factor 2.5. Interestingly, the model predicts a 10 um 
silicate feature with the same intensity as in the observations, 
but with an overestimated continuum. 

That this centro-symmetric disc model does not accurately 
reproduce the observations should not be taken too seriously 
as an argument against backwarmi ng, since the t arget a ppears 
to be a rather complicated object: Malbet et al. (1998) found 
by means of optical long-baseline interferometry that the disc 
should present a hot spot located 15AU from the st ar with a 
flux ratio of 4-5 magnitudes in K; Wan g et al.l (120041) also de- 
tected a companion at 200 AU (0.5") with a similar flux ratio. 
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X (m) k (m) 



Fig. 2. The SED of FUors. Solid line: two-layer disc model including viscous heating and backwarming. Dashed line: black-body 
disc model including viscous heating and backwarming. Dotted line: disc model with viscous heating alone. 



3.2.2. V1057 Cyg 

VI 057 Cyg is well described with an accretion disc featuring 
an accreti on rate M V M = 2.0 X 10~ 5 M|/yr (in the range pre- 
dicted by iKenvon et aljfl988l) . as shown in Fig. |2 (top-right 
panel, dotted line). Backwarming is needed to account for the 
SED in the range 10-100 |xm with a flaring H/r m 0.35 at 
10 AU from the star. The black-body approach either overesti- 
mates the 20-50 um fluxes by a factor 3 (dashed line) or under- 
estimates the 50-100 ^im SED by a similar factor (not shown). 
The two-layer disc model fits the data well in the range 1- 
100 ^im (solid line). It is worth noticing that the model does 
produce a strong silicate feature while the flux in the N-band 
seems to hint to the contrary — it is at the level of the contin- 
uum predicted by the model. 

It is however worth noticing that this disc, though compara- 
ble to FU Ori in terms of accretion properties, presents a much 
larger flaring. If, for some unknown reason, it features a lower 
viscosity, this could account for a larger disc mass and a more 
flared structure. On the other hand, it cannot be ruled out that 
part of the far-IR excess is not to be imputed to backwarming 
and could result from other phenomena, like an envelope or 
self-gravity-related energy dissipation, in which case the flar- 
ing requirement is only an upper limit. 

3.2.3. V1 51 5 Cyg 

V1515 Cyg presents an SED very similar to that of VI 057 Cyg 
at A < 20 |xm but shifted down by a factor 5, which the larger 
distance roughly explains. It is therefore not surprising to find 



a similar accretion rate M+M = 1.3 x 10~ 5 Mg/yr with a stan- 
dard disc model (Fig. [2] dotted line in the bottom-left panel). 
Despite of the lack of data at A > 20 (im, the 10-20 ^m fluxes 
require backwarming with a flaring H/r » 0.35 at 10 AU from 
the star. However, it is not possible to convincingly discrimi- 
nate between the black-body and the two-layer disc models. 

3.2.4. ZCMa 

Z CMa is a close (0.1") binary jKoresko et al J 19891) YSO with 
FUor characteristics. The total flux can be described as emerg- 
ing from an intensively accreting FUor (M ~ 10~ 3 M G /yr) for 
A < 10 \im, yet it fails at reproducing the almost flat far-IR 
SED . A closer look at t he SED of the individual components 
(see (Thiebaut et alJI 19951 and references therein) hint that the 
primary is a cool YSO in its primary phase of evolution (it 
peaks at 5 ^im) while the secondary features an FU Ori-like 
SED. Z CMa B is fairly well described with an accretion disc 
M+M = 1 Ax 10~ 5 M Q /yr, but the lack of data at A > 5 (xm rules 
out checking the backwarming hypothesis (Fig. [5] bottom-right 
panel). 

3.3. Reliability of the model 

This two-layer model proves that backwarming is a credible 
explanation for the mid- and far-IR excess among a few FUors, 
without the need to invoke the contribution of a circumstellar 
envelope. It does not self-consistently tackle the determination 
of the vertical structure, so I discuss here the possible caveats 
linked to these simplifications. 
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Table 2. Disc mass in FUor models: column density S at 1 and 
10AU from the star, cumulated disc mass Ma within 10 and 
100 AU, and radius R„ at which the disc's vertical gravity field 
equals that of the star. Masses are lower values and radii /? s 
upper values, obtained with a high viscosity parameter a = 0. 1 . 





E(10 J kg 


m- 2 ) 


M a 


(M e ) 


R g (AU) 




1AU 10 AU 


10 AU 


100AU 




FUOri 


21 


4 


0.020 


0.54 


29 


VI 057 Cyg 


15 


3 


0.015 


0.32 


50 


V1515Cyg 


10 


2 


0.010 


0.24 


61 


ZCMaB 


11 


2 


0.012 


0.31 


40 



3.3.1 . Simplification of the vertical structure 

Freeing the flaring parameter H/h may seem hazardous, since 
its value strongly affects the amount of irradiation caught 
by the disc and therefore the strength of the IR excess (ex. 
iKenvon & Hartmannlll99l( for backwarmed discs), so a self- 
consistent determinati on of the irradiated surface by more 
elabo rate simulations JP'Alessio et al.lll998t iDullemond et alJ 
2002| may seem more able to corroborate or invalidate the 
backwarming scenario. Yet the latter simulations hide a flar- 
ing hypothesis, since they depend on the up-to-now poorly 
constrained mass column in the disc; in other words, they 
do not help to rule out or support the backwarming scenario 
better than an empirical approach, as far as the diagnosis is 
based on the SED. In the future, they should rely on better 
constrained masses and provide a reliable diagnostic, when 
the Atacama Large Millimetre Array (ALMA) probes the disc 
mass at the lOAU-scale or when IR long-baseline interfer- 
ometry measures flaring-related asymmetries at the AU-scale 
— for instance with AMBER on the Very Large Telescope 
Interferometer (VLTI). 

The two-layer approximat ion may also seem problematic, 
but lDullemond & Natta ( 2003 ) showed that its SED predictions 
stand close to those of full transfer simulation s. On the con- 
trary, it is not fit for spectroscopic studies JAikawa et al.l2 002) 
since spectral lines are very sensitive to the temperature and 
density profiles; in particular, it cannot be used to issue a diag- 
nostic using the strength and shape of the silicate feature. 

3.3.2. Ignoring the disc's own gravity 

FUor disc models predict large masses, as shown in Table [2] 
so that the gravity of the disc can overwhelm the stellar gravity 
in the outer, massive parts of the disc (self-gravitating discs, 
see lBertin & Lodatoll999HHuref200ol and references therein). 
These works single out three effects: 

1. An additional vertical gravitational field that flattens the 
outer parts of the disc and should prevent the m from b eing 
illuminated by the central source of radiation llHurefe OOO). 
which would include backwarming from inner disc regions. 

2. An additional horizontal gravitational field so that the 
disc is no longer in Keplerian rotation. This affects the 
viscous d issipation and the distribution of the matter. 
( Berti n & Lodatol 19991) 



3. Possible instabilities in self-gravitating regions. The disc is 
then more likely to depart from a symmetrical geometry 
and presents an additional dissip ation of en ergy, hence an 
enhanced IR excess. lLodato & Bertinl 12001 ) proposed this 
phenomenon to explain the far-IR excess among FUors. 

PhenomenonHis critical when determining the SED of a back- 
warmed disc. For each model fit, I computed an estimate of 
the radius R g at which the disc undergoes a transition toward 
self-gravity, namely where the vertical gravity field produced 
by the disc at its surface equals that of the central star. As 
shown in Table |2] the three successfully fit FUors should have 
R g ~ 50 AU with the hypothesis of a high viscosity coef- 
ficient a = 0.1. In this case, self-gravity should mostly af- 
fect the coolest regions and therefore the SED at A > 50 \im. 
Phenomena 12131 mostly affect viscous dissipati on and can ac - 
count for the far-IR SED for a ~ 0.01 dLodato & Bertinl200lh . 
It is not impossible that the nature of the IR excess should be 
imputed to contributions of both backwarming and self-gravity. 

It should be emphasised that no study of self-gravity in- 
cluded the influence of irradiation. The latter possesses a con- 
terweight effect: by heating the outer parts of the disc, it in- 
creases the flaring and — at least within the a viscosity pre- 
scription — decreases the amount of matter needed to sustain a 
given accretion rate. Furthermore, it is not clear whether FU Ori 
discs are as massive as expected from standard steady-disc 
models: it might well happen during the outbursts that only the 
innerer parts undergo a sudden increase of the accretion rate, 
while the outer p arts of the disc conserve TTS-like properties 
*Bel1 etalll99.4. 

4. Conclusion 

Using the viscous and irradiated disc model by Lach aume et alJ 
( 2003), it is possible to account for the strength of the mid- and 
far-IR excess among a few FUors using the viscous luminos- 
ity of the inner parts of the disc as an irradiation source. Yet 
the model makes several approximations in the determination 
of the structure: the temperature profile is assimilated to two 
isothermal layers, a surface and an interior, which should not 
prove critical in the SED diagnostic; the flaring of the disc has 
been left as a free parameter, but it proves no less relevant than 
full numerical simulations, that assume the unknown amount 
of material in accretion discs at a scale < 100 AU; more criti- 
cal is my leaving aside the self-gravity of the disc, which could 
"unflare" the disc and prevent irradiation, but the simulation 
of self-gravity in the presence of irradiation has not yet been 
performed, so its influence still remains speculation. 

This work can be seen as a feasibility study for a future 
vertical structure simulation of irradiated discs that I am de- 
velopping us i ng the radiative transfer formalism presented in 
Mal bet et alJ 11200 II) . Such a simulation is needed to issue a 
spectral diagnostic of irradiation (eg. silicate feature in emis- 
sion), that forthcoming IR long-baseline insterferometers will 
ease with their AU-scale resolution for the closest FUors. In 
particular, MIDI on the VLTI will be able to measure the sili- 
cate feature while disentangling the contribution of an envelope 
from that from an irradiated disc, which an SED diagnostic can- 
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not do dVinkovic et alJl2003l) . One also expects to obtain con- 
straints on the flaring with IR closure phases (AMBER on the 
VLTI) and on the disc mass with ALMA; their determination 
is not model-independent and also requires a reliable model. 
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